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ABSTRACT 

The main objective of this paper is the investigation of the singular 
nature of the crack tip stress field in a nonhomogeneous medium having a shear 
modulus with a discontinuous derivative. The problem is considered for the 
simplest possible loading and geometry, namely the antiplane shear loading. of 
two bonded half spaces in which the crack is perpendicular to the interface. 

It is shown that the square-root singularity of the crack tip stress field is 
unaffected by the discontinuity in the derivative of the shear modulus. The 
problem is solved for a finite crack and extensive results are given for the 
stress intensity factors. 


1 . Introduction 

In a nonhomogeneous medium if the elastic moduli are piecewise constant, 
it is known that the stress field around a crack tip terminating at the inter- 
face has a behavior of the form r a where r is the distance from the crack tip 
and -1 <a<0 (see, for example, [1] and [2]). It is also known that if the 
nonhomogeneous medium has elastic moduli which are continuous with continuous 
derivatives the stress state around the crack tips has the standard square- 
root singularity [3], [4]. What has not been studied so far is the effect of 
the discontinuity of the derivatives of elastic moduli on the crack tip stress 
singularity for a crack terminating at the plane of discontinuity. In this 
paper the problem is studied for the simple case of antiplane shear loading 
of an infinite medium in which the shear modulus y is a function of x only, 
dy/dx is discontinuous along the x=0 plane, and the crack lies in the xz plane 
(Fig. 1). The main objective of the paper is to investigate the effect of 
dy/dx on the stress singularity. Hence, it is further assumed that y is an 
exponential function in x, which appears to lead to a relatively simple for- 
mulation of the problem. 
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2. Formulation of the Problem 

The antiplane shear problem for the nonhomogeneous medium shown in Fig. 
1 having the elastic properties 

y (x) = y Q e ex , x>0 ; y(x) = u o e yx , x<0 (1) 

may be formulated as follows: 

3W X 


V 2 w-j + B 9x = 0 , 0<x<«> , 0<y<°° , 

(2) 

3w„ 

v 2 w 2 + y 3x " 0 , -°°<x<0 , 0<y<” 

(3) 

w-j(0,y) = w 2 (0,y) , 

(4) 

°lxz (0 ’ y) = c, 2xz (0,y) ■ 

(5) 

w-j(x,0) = 0 , 0<x<a, b<x<® , 

(6) 

criy Z ( x ,+°) = P(x) , a<x<b , 

(7) 


where it is assumed that the crack surface traction p(x) is the only nonzero 
external load and that because of symmetry it is sufficient to consider the 
problem for y> 0 only. Expressing the solution of (2) and (3) as 

00 00 

1 " * o ’ 

w l( x,y ) = 2^ f-|(y,a)e _1aX da + - g 1 (x,a)sinyada , (8) 

J ■ J 


oo 

2 

w 2 (x,y) = — g 2 {x,a)sinyada 


we obtain 


f-j = A(a)e my , g-j = B(a)e nx , g 2 = C{a)e A) 
m = - Va 2 +i 8a , n = -a-j - , X = a 2 - » 



(9) 
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The stress components are given by 


a lxz = V 


a 2xz = y o e 


yx 


3W-| 


e 8X 8W 1 

3X 

5 a lyz = 

y o e 3y • 

3W 2 


„ P Y* 2 

3X 

c 

ro 

*< 

N 

y o e 3y ’ 


( 11 ) 


( 12 ) 


Substituting from ( 8) - ( 1 2 ) into the continuity conditions (4) and (5) it 
may easily be shown that 


1 


C(ct)-B(a) = 


[ aA(p)dp 

a^+p^+i Bp . ’ 


(13) 


nB(a)-XC(a) - i [ . 


04) 


Defining now 
g(x) = w ] (x,0) , 


(15) 


from (8), (10) and (6) we find 
b 

iaA(a) = -[ g(t)e l0lt dt , 


(16) 


g ( t ) dt = 0 . 


07) 


By substituting from (16) into (13) and (14) and by using the residue theorem 
to evaluate the improper integrals and the condition (17), the unknown functions 
B and C may be obtained as follows: 


f ^ a(a-j -<*2 + ^“) 

B(a) - ^Cx-nH^-e/^) J 


g(t)e 


t(a r e/2) 


dt , 


(18) 
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C(ct) 


a 

(x-n) (a-|-e/2) 


f u -t(a,-e/2) 
g(t)e 

a 


dt 


(19) 


Thus, g(t) is the only unknown in the problem which may be determined from 
the remaining boundary condition (7). From (7), (8), (10), (11), (16) and 
(18) it then follows that 


o lyz ( x >°) = P(x) = u 0 


e ex 1 


b 

* 

[k^(x,t) + k 2 (x,t)]g(t)dt , a<x<b , 
a 


( 20 ) 


k,(x,t) = lim 
1 y-H-0 


i_ 

2 


rn(g) 

a 


e my e ia(t-x) da f 


( 21 ) 


k ? (x,t) = lim e (t~ x )&/2 
y-++0 


■ a 2 (a-t-a 2 + ) _ (t+x)a, 

— g— e cosay dx . 

J Q a-| (x-n) (a-j - |-) 


( 22 ) 


Referring to (10) and the regularity conditions at infinity it may be noted 
that Re(m)<0, Re(n)<0 and Re(x)>0. The singular behavior of the kernels k.j 
and k 2 may be obtained from the asymptotic analysis of the integrals in (21) 
and (22). By observing that for large values of a, m(a)->- |a | , from (21) we 
find 


k,(x,t) = lim 
1 y-»+0 


2 


co 


00 


e‘l“l y e 1a ( t_x) d a 


00 



— 00 


(23) 


where because of uniform convergence, in the second integral the limit has 
been put under the integral. By evaluating the first integral k-j i s obtained 
as follows: 

k-j (x,t) = + h-j (x,t) , (24) 
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oo 


( 25 ) 


h,(x,t) = Im{[ ('/T+TI - De^^^da} . 

I J a 

o 

where the function h-j is bounded for all values of x and t in the closed 
interval [a,b], a_>0. 

Similar asymptotic analysis would show that the kernel k2 is bounded 
for a>0 (for which t+x>0), and has only a logarithmic end point singularity 
for a=0 (i.e., for x+t->0). We first write 


k 2 (x,t) 


= e 


( t-x) 6/2 


l<2(x,t,a)cOSciy da 


K 2 ( X jt ,a ) 


a 2 (c.i-« 2 tlji) (t+x)a 
ct i (x-n)(a-|-6/2) 


(26) 


(27) 


The limiting behavior of the integrand K 2 for a-*0 and for a-*» is 

K 2 K 2o (x,t,a) = a 2 e“ (t+x)6/2 , for a+O , (28) 

K 2 + K 2 Jx,t,a) = 1 e " a(t+x ) , for «■*» . (29) 


We now express 

k 9 (x,t) = lim 
L y->+0 


k 2 as 


e (t-x)B/2j-f 

o 


K 2o (x,t,a) COSay da + 


+ 


00 

K 2co (x,t,a)C0Sayda +6] 
N 


N 

l< 2 ( x ,t ,a ) COSay da 

E 


(30) 


where e is a very small and N is a very large constant and the constant 6 may 
be made as small as we please by selecting e sufficiently small and N suffi- 
ciently large. Note that in E£a_<N K 2 is bounded and hence, the second integral 
in (30) is finite for all x>0, t>0. Equation (28) shows that the same is true 
also for the first integral in (30). The third integral is the exponential 
integral which may be expanded as 
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r -a( t+x) 

lim cosayda = lim Re{Ei[(t+x-iy)N]} = -log(t+x)-logN-Y 

y++0 j N a y->0+ 

N 2 , .2 m3 3 , 

- 272y (t+x) +3V (t+x) - ••• (31) 

where y is the Euler's constant. One may note that if N is nonzero and finite 
' o 

then the exponential integral in (31) and consequently the kernel k 2 is bounded 
for all x and t in 0<(x,t)<b and has only a logarithmic singularity at x=0=t. 
Since the kernel log(t+x) is square integrable, it may thus be concluded that 
the dominant part of the integral equation (20) has only a simple Cauchy ker- 
nel for a=0 as well as a>0, and the solution is of the form [5] 


g(t) = 


G(t) 

/(t-a) (b-t) 


(0 £ a < t < b) 


(32) 


3. Stress Intensity Factors 

For a>0 we observe that (20) gives °-jy 2 ( x >0) for 0<x<°° and the function 
F defined by 

b 

* 

F(x) = [h-j (x,t) + k 2 (x,t)]g(t)dt (33) 


is bounded in the closed interval a<x<b. Thus (20) may be expressed as 

b 


h e ' BX ^ (x > 0) • I 


- -Ktjdt + F(x) s ( 0<x<co ) 
a (t-x)X + (t) 


(34) 


where 


X(z) = /(z-a)(z-b) , X + (t) = -X"(t) = /(t-a) (t-b) , (z=x+iy) 


(35) 


Consider now the function 


if^rfey 

C 


(36) 
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where the contour C encircles the crack and z is outside C. By shrinking 
the contour to the cut, from (34) -(36) it may be shown that 

e" 8x a-jy Z (x,0) = - 4 ) (x) + F(x), (0<x<a, b<x<~) . (37) 

On the other hand, from (36) it follows that [5] 

4>(z) = f[fj " P(z) (38) 

where P(z) is the principal part of G/X at |z| = ». From (37) and (38) we 
find 

7" e " ex a iy Z (x ’ 0) = _ Iffy + P(x) + F(x)} (0<x<a ’ b<x<ro) • (39) 

By observing that 

X(x) = /(x-b)(x-a) , (x>b); X(x) = -/(a-x)(b-x) , (x<a) , (40) 

the Mode III stress intensity factors at the crack tips may now be defined 
and expressed as follows: 

k,(b) = 11m v^Tx^bl a,...(x,0) = -y_e eb — G - ^ - 
3 x+b lyz 0 /(b-a)/2 

= -11m y(x)/2(b-x)g(x) , (41) 

x->b 


ko(a) = 11m /2(a-x) a, (x,0) = y e Ba G ^ — 

3 x->a lyz 0 /[b-iT/2 

= lim y(x)/2(x-a)g(x). (42) 

x->a 

For a=0 we define the stress intensity factor at the crack tip x=0 as 
follows: 

k-(0) = lim /^2x a 9 (x,0) . (43) 

J x-»-0 zyz 
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To calculate k 3 (0) the asymptotic analysis of around the crack tip 
x=0 is needed. From (9), (12) and (19) cr 2 y Z ma y & e expressed as 


o 2yz (x,y) = U() e* x § 


X x 

C (a ) e aCOSayda 


o 

b 


xX l f h,(x,t)g(t)e (Bt - xx)/2 dt 


U.6 I “o 

0 TT | 3 


(44) 


where 


h^(x,t) = lim 2 
J y-^+o 


a 0 X-a-, t 
2- 1 


a^e 


U- n)( or e/2) y *• < 45 > 

By observing that a-j->a, a^-m, x->a and n->-a for ct-*», from (45) we obtain 


h~(x,t) = lim 
y ^+0 


r -a( t-x) , 
e ' COSay da 


a 0 X-a n t 


r2a 2 e „-a(t-x)n j 

C (x-n)(a r e/2) ■ e ]da 


(46) 


where the second integral is uniformly convergent. Note that in (46) t>0, x<0 
and t-x>0. Thus, from 


lim 

y-vt-0 




y-*+u 


(47) 


and (46), a 2yz may be obtained as follows: 


a 2y Z ( X *®) ~ 


yx/2 1 


u 

' H 3 (x,t)]g(t)e 8t/2 dt , 


(48) 


where the bounded kernel is given by the second integral in (46). 

If we now substitute from (32) into (48) and follow an analysis similar 
to (33)- (40), the asymptotic expression for a 2 y Z ma y be written as 


a o .„(x,0) = y e yx/2 G Cx}e e . X/ + K ( x) 

2yz 0 /r-x)(b-x) 


(49) 
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where K is a bounded function. Thus, from (49) and the definition of the 
stress intensity factor at x=0 as given by (43) it follows that 

k-(0) = u = lim w (x)<Sx g(x) . (50) 

d 0 ^72 x->0 

The results obtained in this section clearly show that the square-root 
character of the crack tip singularity is unaffected by the discontinuity in 
the derivative of the shear modulus. 

From (10) and (22) it may be seen that for y = S (i.e., for the case of 
single nonhomogeneous plane) the kernel k 2 is identically zero and the expres- 
sions (24) and (25) defining k-j remain unchanged. Thus, with k 2 =0, (20), 

(24) and (25) give the density function g(x) as defined by (32) and the stress 
intensity factors are then obtained from (41) and (42). 


4. The Rigid Half Space 


In the special case of an elastic half space bonded to a rigid half 
space, w 2 =0, (8) is still valid and A and B are given by (16) and (13) (with 
C=0), respectively. Following an analysis similar to that of Section 2 (equa- 
tions ( 20 ) — ( 31 ) ) , it may easily be shown that in this case the integral equa- 
tion (20) becomes 

b 


P(x) = y Q e 


Bx 1 

TT 


¥ x) - 


J_ -B(t-x)/2 

t+X 


a 


- h 2 (x,t)e^ t-x ^ 2 ]g(t)dt , (a<x<b) 


(51) 


where h-j is given by (25) and 


h 2 (x,t) 


00 



a l (“l 

0 1 1 



-C*1 (t+x) 

e 


e ~a(t+x) 


]da . 


(52) 


We note that for B=0 h^ and h 2 vanish and (51) would reduce to the following 
known integral equation for the homogeneous half plane x>0 for which w(0,y)=0: 


— p(x) = 1 
y o 71 


‘t-x 


^)g(t)dt , (a<x<b) 


(53) 
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For a>0, clearly the solution of (51) is of the form (32) and the solu- 
tion may easily be obtained by following the technique described in, for example, 
[6]. In the limiting case of a=0, by expressing the solution of (51) as 

fit “i 

e y g(t) = G(t) (b-t) 't , (-l<Re(co-, ,aj 2 )<p) (54) 

and by following the function theoretic method, it may easily be shown that 
(see, for example, [5] or [6]) 

G(b)cotiTaj-j = 0 , (55) 

G(0)(cotuaj 2 - $1 - n7I ^) = 0 . (56) 

Equations (55) and (56) are identical to the characteristic equations which 
would result from (53) and give = -1/2 and = 0 as the acceptable roots. 

It is, therefore, seen that at t=0 g(t) is bounded. Physically, this result 
is indeed expected, as the half plane problem with w-,(0,y) = 0 corresponds to 
the antisymmetric problem for the infinite plane for which y(x) = y Q e 1 1 and 
which has a crack on the x axis along -b<x<b subjected to antisymmetric shear 
tractions o yz (x,0) = p(x) = -p(-x), (-b<x<b). 

5. Results and Discussion 

The calculated results are shown in Figures 2-4 and Tables 1 and 2. 

Figure 2 shows the stress intensity factors for an infinite plane with a 
shear modulus y(x) = y Q e subjected to uniform shear tractions a-| yz (x,0) = 

-p . These results are analogous to those given in [4] for the pressurized 
crack. As in [4], the normalized stress intensity factors are independent 
of the crack location d (see Fig, 1) and are functions of the dimensionless 
parameter gc only. At first sight the result given in Fig. 2 to the effect 
that the stress intensity factor k^(a) onthestiffer side of the medium is 
greater than k^(b) may appear to be somewhat paradoxical. However, this result 
may easily be explained by considering the corresponding crack surface dis- 
placements given in Fig, 3 and obtained from 
x 

w-|(x,+0) = j g(t)dt , (0<x<b) . (57) 

o 
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The figure shows the normalized crack surface displacement w-j(x,+0) for the 
nonhomogeneous infinite medium with y(x) = y Q exp(ex), ec=-2, and for the 
homogeneous planes with shear moduli y=y(0)=y Q , y=y(c)=y Q exp(-2) and y=y(b)= 
y Q exp(-4). In the homogeneous planes the stress intensity factors are indepen- 
dent of y and are given by 


k 3 (0) = k 3 (b) = p Q /c 


(58) 


whereas the crack surface displacement is inversely proportional to y, i.e. 


w(x,+0) = /x(b-x) , (0<x<b) 


(59) 


From Fig. 3 it is seen that near the crack tip x=0 the crack surface displace- 
ment w(x,0) for the nonhomogeneous medium is greater than that for the homo- 
geneous medium having the modulus y(0) = y . Since the stress intensity fac- 
tor k 3 is related to the magnitude of the crack surface displacement derivative 
(see (41), (42) and (50)), it would, therefore, be expected that (for e<0) at 
x=0 k 3 for the nonhomogeneous medium would be greater than p Q /c, the value 
for the corresponding homogeneous medium. Even though near x=b w(x,0) is 
considerably greater than that near x=0, it is still smaller than the displace- 
ment for the homogeneous medium having y=y(b) which also has k 3 (b) = p Q /c. 

With (41), this would then explain the trend for the stress intensity factor 
at x=b shown in Fig. 2, namely that k 3 (b)<p Q /c. 

Some sample results for the crack surface displacement in bonded nonhomo- 
geneous half planes with a crack along 0<x<b are shown in Fig. 4. 

For uniform crack surface traction a iy Z ( x >0) = “P o Tables 1 and 2 show 
the calculated stress intensity factors normalized as 


k 3 (a > 

P 0 /f 


k(a) 


k 3 (b) 

P 0 ^ 


= k(b) 


(60) 


Aside from the validity of the general trends given in Fig. 2 for the infinite 
medium, these results also conform to the broad principle that in bonded non- 
homogeneous solids as the crack tip approaches the interface the corresponding 
stress intensity factor tends to decrease if the crack lies in the medium 
with the smaller shear modulus and increase if the crack lies in the stiffer 
medium. 
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It should perhaps be pointed out that even though the main result of 
this study namely, that the square-root nature of the stress singularity is 
unaffected by the "kink" in u(x), is based on a specific choice of the shear 
modulus (y Q e^ x for x>0, y Q e yX for x<0, clearly the conclusion should be 

valid for any continuous function y(x). One may, therefore, also conclude 
that the expression 

MO) = lim y(x)v^x — w(x,0) 

3 x-sO 8X 

for the stress intensity factor at a crack tip x=0 is valid for any y(x) which 
is continuous at x=0. 
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Table 1. Stress intensity factor ratio k 3 /p 0 v^ in bonded 
nonhomogeneous half planes. 
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Table 2, Stress intensity factor ratio kg/pg/c in 
bonded nonhomogeneous half planes. 
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Fig. 1 


Geometry for bonded nonhomogeneous half planes. 



Fig. 2 Stress intensity factors for an infinite nonhomogeneous plane 
subjected to uniform crack surface traction a^ z (x,0) - -p Q . 



Fig. 3 Crack surface displacement in an infinite nonhomogeneous plane 
under uniform crack surface shear loading c (x,0) = -p Q ; 
shear modulus y(x)=y 0 e ex , gc=-2. y 


y 



x/c 


Fig. 4 


Crack surface displacement in bonded nonhomogeneous half 
planes under uniform antiplane shear loading a-jy Z (x,0) = 
-p 0 ; u(x)=u 0 exp(ex) , x>0; y(x)=y 0 ex p(yx), x<0, B=-l. 
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